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Abstract. The spectroscopy of electron capture neutri- 
nos emitted from nearby pre-supernova collapsing stars 
before the neutrino trapping sets in, can yield useful in- 
formation on the physical conditions and on the nuclear 
composition of the core. The neutrino spectrum depends 
on the thermodynamic conditions of the core, the nu- 
clear abundances, the lepton fractions and relevant nuclear 
properties. In the pre-trapping core, the density ranges 
from 0.1 — 100 10^" g/cm'^ and the temperature from 
0.2 — 1.5 MeV. The nuclear abundances can be obtained 
under the assumption of Nuclear Statistical Equilibrium 
(NSE). The nuclear abundances as well as the electron 
capture rates are thus determined, among other things, by 
the nuclear binding energies and the free nucleon chem- 
ical potentials. Because shell and pairing effects persist 
strongly up to temperatures of ~ 0.5 MeV, any equa- 
tion of state (EOS) relevant to this phase of the collapse 
must reproduce well the zero temperature nuclear proper- 
ties and it must show a smooth transition to the known 
high temperature and high density limits. In this work we 
use the microscopic Relativistic Mean Field (RMF) theory 
based on a Lagrangian with non-linear self-interactions of 
the (T-field for the neutron-rich nuclei of interest in the 
f — p shell to determine nuclear chemical potentials. We 
compare these results with those computed from an EOS 
calculated with the macroscopic liquid drop model. We 
also discuss extensions to finite temperature and we in- 
corporate nuclear lattice effects into the microscopic cal- 
culations. 

Key words: Stars: supernovae: general - Equation of 
State - Nuclear reactions, nucleosynthesis, abundances 



1. Introduction 

The collapse of the core of a massive star is believed to 
lead to the explosion of a type II supernova. So far, how- 
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ever, numerical simulations of such a collapse have not 
been able to demonstrate convincingly the explosions ob- 
served in nature. Such simulations require many detailed 
and complex physical inputs, such as the equation of state 
of dense and warm stellar matter, the hydrodynamics of 
the shock wave, or the energy transport by neutrinos. The 
calculation of the equation of state of the stellar core ma- 
terial at densities below the nuclear saturation densities 
(10^ < p < 10^^ g/cm'^) constitutes, technically, one of 
the most computer intensive parts, because it requires a 
high accuracy due to its influence on the early part of 
the supernova development (see e.g. Cooperstein & Baron 
1990). 

Initially, at densities of 10^ g/axi^ the core is com- 
posed of heavy nuclei immersed in a electrically neutral 
plasma of electrons, with a small fraction of drip neutrons 
and an even smaller fraction of drip protons. Nevertheless, 
the overall density being much lower than nuclear matter 
density, the average volume available to a single nucleus is 
much greater than that of the nuclear volume. Under such 
circumstances the equation of state can be approximated 
by that of a Boltzmann gas of heavy nuclei, and drip par- 
ticles in Nuclear Statistical Equilibrium (NSE), together 
with a highly degenerate, electrically neutralizing gas of 
non-interacting electrons. This state of affairs persists up 
to densities of ^ 10^^ g/cm'^, where the ratio of the 
nuclear radius {Rn ~ A^^^) to the Coulomb interaction 
radius {Re - {A/p)^/^) given by Rn/Rc ^ lO^p^/^ is no 
longer small ^ 1 and the nuclei start interacting with 
each other by means of Coulomb lattice effects. Finally, 
at values around 10^^ g/cm'^ the density of the "vapor" of 
free nucleons surrounding the droplets of nuclear matter 
is too high to justify the treatment of nuclei as a Boltz- 
mann gas of non-interacting particles (see e.g. Hillebrandt 
(1994) and references therein). 

The equation of state of moderately neutron-rich mat- 
ter at intermediate temperature and density can substan- 
tially influence the energy spectrum and the flux of the 
early phase neutrinos from electron capture on nuclei and 
free protons. Should the neutrinos from the collapse phase 
be detected from a nearby supernova in an underground 
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neutrino detection experiment such as the Sudbury Neu- 
trino Observatory (SNO) or the Super-Kamioka (SK), 
such spectra could be possibly measured. This will be pos- 
sible only under the assumptions, that such a supernova 
explosion will occur at a distance of approximately 1 kpc 
and that these new generation of neutrino detectors with 
increased sensitivity are available at that time (see Su- 
taria & Ray (1997), hereafter SR). The neutrinos from 
the early phase of the collapse, before they are trapped by 
dense overlying matter are especially interesting for sev- 
eral reasons. They bear information directly related to the 
conditions under which they are produced (SR). The pre- 
trapping phase of the emission of core-neutrinos continues 
until the overlying density grows to about 10^^ g/cm'^. 
The emitted neutrino spectrum depends upon the indi- 
vidual Q- values Qi of the nuclei involved in the electron- 
capture reaction. Both these Q- values as well as the overall 
abundance of particular nuclei depend upon the binding 
energies of these nuclei. As the electron capture on pro- 
tons has different neutrino emission characteristics than 
the capture on nuclei, the emitted spectra will be differ- 
ent for different admixtures of these components. The nu- 
clear equation of state (EOS) controls the heights and the 
centroids of the neutrino spectrum such as the example 
shown in Fig ol This example was computed using the 
EOS, in particular the expressions for nuclconic chemical 
potentials as given by Fuller (1982). 

During the collapse, it depends also on the entropy of 
the core whether most of the electron captures occur on 
free protons or on heavy nuclei. At relatively low entropies 
(e.g. S/Ub ~ 1) captures on heavy nuclei dominate the to- 
tal rate. Such entropies of the core do occur for stars of 
main sequence mass 10 - 25 Mq . As the entropy connected 
with nuclear excitations and the entropy of the free nu- 
cleons form a significant part of the total entropy of the 
stellar core, this forms another reason for evaluating the 
nuclear properties of the stellar core with accuracy. The 
relevant temperature range of the stellar core is approxi- 
mately 0.2 to 1.5 McV at such low entropies for a density 
range of 10^ - 10^^ g/cm^. 

Investigations on the equation of state of supernova 
matter have by now spanned several decades (Bethe et 
al. 1979, hereafter BBAL; Lamb et al. 1978, hereafter 
LLPR78; Lattimer et al. 1985, hereafter LPRL85; Lamb 
et al. 1983, hereafter LLPR83). At the same time, ever 
since the discovery of neutron stars as pulsars more than 
thirty years ago, work on the equation of state of "cold, 
catalyzed" matter, i.e. matter at zero temperature and 
under thermodynamic equilibrium, has been pursued by 
a number of scientists (see Pethick & Ravenhall 1995 and 
references therein). The equation of state of hot super- 
nova matter has been in many ways an extension of the 
early work done on the EOS of neutron star matter at zero 
temperature (Baym et al. 1971a (BPS); Baym et al. 1971b 
(BBP)). Since these investigations were based on the liq- 
uid drop mass formulae of Myers & Swiatecki (1966), the 



calculations of the equations of state for supernova have 
also followed largely this semi-empirical approach. More 
recent attempts of the EOS of supernova matter by Swesty 
et al. (1994) are based on the finite temperature, com- 
pressible liquid droplet model which attempts to take into 
account the nucleon-nucleon interaction in analogy to the 
density dependent Skyrme type energy functionals. On 
the other hand, there are microscopic investigations based 
on the Hartree-Fock or on the Hartree-Fock-Bogoliubov 
scheme with effective nucleon-nucleon interactions, such 
as that of Negele & Vautherin (1973) above the neutron 
drip point or that of Haensel et al. (1989) below the neu- 
tron drip point. 

However, as far as pre-supernova conditions are con- 
cerned, the liquid-drop (or droplet) based versions of the 
EOS cannot be applied so easily because in this temper- 
ature range (0.2 to 1.5 MeV) shell and pairing effects 
dominate the nuclear structure. It is only at tempera- 
tures beyond ~ 1.0 to 1.5 MeV, that nuclei can be dealt 
within the liquid drop approximation. Thus an equation 
of state which is used to calculate nuclear properties in 
the pre-supernova stage should take into account "micro- 
scopic" shell and pairing effects. At higher temperatures 
and densities, well after neutrino trapping at 10^^ g/cm'^ 
they should smoothly go over into the high-temperature 
high-density EOS based on the liquid drop model of nu- 
clear matter. Furthermore the computation of the EOS by 
theoretical methods needs to be supplemented by experi- 
mental data and their trends. 

In this context, we note the work of El Eid & Hille- 
brandt (1980) (hereafter EH) which computed the equa- 
tion of state of supernova matter and made use of the 
so called droplet model (DM) and its extension. The DM 
model was an early attempt which went beyond the con- 
ventional liquid drop model (LDM) also due to Myers 
& Swiatecki (1966) with the inclusion of a shell correc- 
tion term due to von Groote et al. (1976). The LDM was 
used for the EOS computations of much of the subsequent 
works mentioned below. The Droplet Model is an analyti- 
cal extrapolation of the Thomas Fermi model. Von Groote 
et al. also gave a new fit to the DM parameters to all nu- 
clear masses known at that time. In this paper, we present 
a microscopic calculation based on a suitably generalized 
Relativistic Mean Field theory adapted to the stellar col- 
lapse problem (we refer also to Sutaria (1997), Ray et al. 
(1999) and Sutaria et al. (1997)). This automatically ac- 
counts for the nucleon-nucleon interactions and for shell 
and pairing effects, and as we show later reproduces the 
currently known nuclear masses with substantial accuracy. 
A microscopic approach like this is expected to be more 
reliable for nuclei relevant in the astrophysical context, 
which are often very neutron rich and are thus far from 
the valley of /^-stability. 

In section we summarise the existing equations of 
state based on the liquid drop approach; this helps to es- 
tablish the notation and to compare with previous ap- 
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proaches. In section ^ we describe the Relativistic Mean 
Field theory calculations as adapted to the supernova con- 
text. In this section we also present the results of these 
investigations. In section H we compare our results with 
those of liquid drop based models and compare with the 
earlier work of Cooperstein and Baron (1990) and refer- 
ences therein. In section ra we describe the generalization 
to finite temperature and discuss the modification of our 
results due to low but finite temperature effects. In the 
concluding section |g we make a comparison with the re- 
sults of EH and ours at the initial stage of stellar collapse 
and summarize our conclusions. 



2. Suitable Equations of State for each Stage of 
the Collapse 

The initial efforts to compute the EOS, and indeed most 
of the literature available so far on the EOS of super- 
novae concentrate mainly on matter at densities beyond 
10^'^ g/cm'^. Most of them are based on the liquid drop 
model of the nucleus, as the initial work of BBAL, LLPR 
and their co-workers. Later finite temperature extensions 
of the existing EOS as for instance the work of Baron et 
al. (1985) (hereafter BCK) and Cooperstein (1985) took 
into account the nuclear compressibility. The most sophis- 
ticated attempts in developing an equation of state at 
higher temperatures and higher densities, close to nuclear 
saturation densities, have been finite temperature Hartree 
Fock calculations, as those of Bonche & Vautherin (1981). 
These Hartree-Fock approaches rely on the spherical ap- 
proximation and on phenomenological models for the nu- 
clear potentials, usually of Skyrme type, which are fitted 
carefully to properties of finite nuclei. 

As discussed in Sect. || the composition of stellar mat- 
ter at densities between 10^ and 10^^ g/cvo? and at tem- 
peratures between 0.2 and 1 MeV can be approximated by 
a 'soup' of nuclei including drip neutrons, protons, alpha- 
particles and an ensemble of heavy nuclei. The nuclear 
abundances in this pre-supernova stage with a density 
p « lO'^ g/cm^ peak in two regions of the isotope table, 
one around the silicon group and the other around the iron 
group. The nuclear species in each of these groups remain 
in quasi-equilibrium among themselves but the groups of 
nuclei are not in quasi-equilibrium with each other at the 
lower densities (Hix & Thielemann 1996). The fraction of 
drip nucleons remains low during the collapse until nuclei 
merge into a uniform nucleon sea, possibly through a bub- 
ble phase. This is because the entropy per nucleon of the 
system remains close to 1 (in units of k^) implying a high 
degree of order in the system (BBAL). 

In the macroscopic incompressible liquid drop ap- 
proach to the cold EOS above the neutron drip (BBAL 
1979; BBP 1971), the ensemble of heavy nuclei and drip 
nucleons is considered on a lattice of number density njv 
consisting of a single heavy nucleus with mass and charge 
numbers (A,Z) having a nuclear volume Vm immersed in a 



sea of drip neutrons of density n„ and electrons of density 
rig. For such a system, the effective free energy density 
Etot at zero temperatures (the internal energy density of 
the system), is given by BBP: 

Etot[A,Z,nN,VN,nn) = 

riNiWN + Wl) + (1 - VNnN)En{nn) + E,{ne). (1) 

Wl is the Coulomb lattice energy of the system, (1 — 
VNn^) is the fraction of the total volume occupied by 
the neutron gas. En is the energy density of the non- 
relativistic free nucleon gas consisting primarily of drip 
neutrons and Ee{ne) is the energy density of the relativis- 
tic electrons. 

The lattice energy of the system, assuming a bcc lattice 
structure, is expressed as Wl = — 1.82Z^e^/a, where a = 
(2/nAr)^/'^ is the lattice constant. The quantity Wn is the 
nuclear matter energy of a nucleus in equilibrium with the 
free neutron gas, and is given by BBP: 

Wn = [{l-x)jnn+xmp]c^A+[W + WtMck + We,^ch\A 



W, 



surf 



A3 + Wc.oZ^A- 



(2) 



where x is the proton fraction of the heavy nucleus, m„ 
and TOp are neutron and proton masses, W is the energy 
of bulk nuclear matter per nucleon, Wsurf is the coeffi- 
cient of the surface energy term and Wcfi is the nuclear 
Coulomb energy, which is taken to be the Coulomb energy 
of a uniformly charged sphere with charge Z. There are 
two corrections to the total Coulomb energy; a term due to 
the finite surface thickness which is denoted by WthickA, 
and another due to the proton exchange term WexchA, 
but both these are generally negligible (BBP). The total 
nuclear energy Wjv, including the corrections due to lat- 
tice energy term (because the system consists of a lattice 
of positively charged nuclei sitting in a sea of electrons) 
and the 2nd order corrections to this lattice energy (due 
to the non-zero Coulomb radius) is given in terms of of the 
proton fraction x, the density in the nuclear interior pjsi, 
the nuclear volume Vn and the packing fraction i.e. the 
fraction of total volume occupied by nuclei {u — p/ pn), 
as (Bethe et al. 1983, hereafter BBCW): 



Wn{x,pn-,Vn-,u) 



Wbulk 
„2„2 



2QQx^{l-xfA^ 



+Px'pi,VS{l-^u^ + ]^u), 



(3) 



where the factor j3 ~ 3/5, and Wbuik is the energy of 
infinite nuclear matter without surface or Coulomb, but 
with asymmetry effects, given in BBCW: 



WbuL 



[-16 + 29.3(1 -2a;)2], 



(4) 



The above expression for Wbuik is based on the simple 
incompressible liquid drop mass formula. The first term is 
the volume energy of the nucleus per nucleon. The second 
term represents the nuclear volume asymmetry energy and 
the coefficient of the volume asymmetry energy S^ . Here 
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we use the value Sy = 29.3 obtained in the literature from 
fits to the experimentally known nuclear masses (Swesty 
et al. 1994). 

The mass number A of the most dominant nuclear 
species is determined by a minimum of the free energy 
with respect to A at constant x, unA, tinVn, rin. This 
leads to the conditions that the most stable nuclear system 
(for a given e~-fraction Ye = 2;, drip neutron density n„, 
packing fraction u and heavy nuclei density p^r ~ tlnA) 
is that for which the surface energy is twice the Coulomb 
energy: Wsurf = 2WcouiA. 

The chemical potential of the neutrons /i„ is defined 
as the change in the nuclear energy when a single neutron 
is added to the nucleus, keeping the number of protons 
constant. Therefore /z„ can be obtained from Eq. (0) by 



1 



fJ-n 



dW, 



N 



x_ f dWN 
Vn \dpN J,^^^^ A\ dx 



(5) 



Pn ,Vn .v. 



and similarly, the difference of neutron and proton chem- 
ical potentials is obtained by: 



11 = Hr, 



Mp 



1 fdW, 



N 



A\ dx 



(6) 



PJV.VJv," 

Using Eqs. (||) and (^ along with Eq. (||) and the condi- 
tion that the surface energy is twice the Coulomb energy 
for the "mean nucleus", one can obtain relations for //„ 
and fl. Fuller (1982), quoting a private communication 
from Lattimer (1980) gives an expression for fin as :- 



fi„ = -16 + 125(0.5-3;)- 125(0.5 -a;)^ 
, xHl~x)H3^7x) 
A^/^2{l-x) 

and for fi as (BBAL) :- 

290 



(7) 



A = 250(0.5 -x)- 1^.2(1 - 2xr (^ + ^\ 



.,1 , 21-2.^^^^^ 



Note however the comments pertaining to these, dis- 
cussed in the next section. This model described in Eqs. 
(I - |) (cited as BBAL/BBP-EOS later on) is based on 
an incompressible liquid drop. It uses the surface energy 
term of Ravenhall et al. (1972) and it assumes that the 
actual nuclear density can be approximated by the nu- 
clear density at saturation. The relations for /i„ are taken 
from Fuller (1982) and are meant for near-symmetric nu- 
clear matter for which x ~ 0.5. The correction terms in 
fin are discussed in the next section. In Fig.^ we display 
the LDM expression of Eq. u\ for the chemical potential of 
the neutrons in the chain of Mn-isotopes along with the 
values obtained form relativistic mean field calculations 
Mn|_RMF which will be discussed later. Similar diagrams 
for the elements Fe, Co, Ni, Cu, Zn, Ga, and Ge can 
be found in Sutaria (1997). For the curve labeled BCK 
in Fig.g, modifications were made to take into account 
the difference between the actual nuclear density and the 



saturation density in the compressible EOS developed by 
Baron et al. (1985) (BCK) and Cooperstcin (1985) (cited 
as BCK-EOS lateron). 

The BCK-EOS differs from the BBAL/BBP-EOS in 
three points: (a) a term is incorporated in the expres- 
sion for Wn/A which takes into account the nuclear com- 
pressibility, (b) the actual nuclear density po may differ 
from the saturation nuclear density of infinite nuclear 
matter ps, (c) pg is allowed to depend on the proton 
fraction x, where the functional form of this dependence 
{(t>{x) = p,{x)/0A6 fm-3 = 1 - 3(0.5 - x)^) was obtained 
by fits to the density dependent HE calculations with the 
force SKM of Bonche & Vautherin (1981). In the com- 
pressible model, the expression for Wn has the form: 



W, 



N 



-16 



1 



29.3 (l-2x)^ -1- —K{1 



+ 75.4x^(1 -a;)30~3( 



18 



1 



-U3 



(9) 



Where K is the coefficient of nuclear incompressibility, 
9 = pa/ps is the ratio of the actual nuclear density po 
to the saturation nuclear density Ps, and u is the pack- 
ing fraction of the nuclei, such that u = p/ Pq. The term 
dependent on the incompressibility was added in order 
to obtain an accurate description of the phase transition 
into the bubble region without producing any discontinu- 
ities in the pressure or the temperature variables in the 
numerical calculations simulating the core collapse. The 
low entropy liquid drop-based equations of state (hereby 
the LDM-EOS) all apply to warm dense matter. The ap- 
proaches vary: BBAL use the phenomenological approach, 
which included many basic features of the microscopic cal- 
culations of LPRL. BCK use the "top down" route, where 
the relevant nuclear parameters like bulk incompressibil- 
ity, saturation density including its dependence on Yg, on 
the nucleon effective mass, or the bulk asymmetry energy 
coefficients, etc. were taken as input quantities for the cal- 
culation of thermodynamic variables, rather than as out- 
put for a specific set of nuclear force parameters, as in 
the microscopic calculation. The BCK series of EOS were 
constructed for numerical applications in supernova hy- 
drodynamics and therefore with physical and numerical 
simplicity as important considerations. The composition 
of warm dense matter is determined under the conditions 
of mass and charge conservation, by the pressure balance 
at the surface of a heavy nucleus surrounded by light nu- 
clei like a-particles and by a " vapour" of free protons and 
neutrons, and under the assumption that the strong and 
electro-magnetic reactions have reached nuclear statistical 
equilibrium (NSE) (Cooperstein & Baron 1990). 

An explicit dependence on temperature was intro- 
duced in the BCK-EOS free energy (F) term by using 
the expression: 



F = Wn- -t T 

A m 



(10) 
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where m* /m is the effective mass of the nucleons in the 
nucleus and a is level density parameter in the Fermi-gas 
model of the nuclei. 

3. The Relativistic Mean Field (RMF) Approach 

The Hartree-Fock (HF) or Hartree-Fock-Bogoliubov 
(HFB) methods of calculating nuclear properties depend 
on the knowledge of a suitable nuclear potential. Usually 
Skyrmc or Skyrme-type forces are used for this purpose. 
In recent years the microscopic description of ground state 
properties of finite nuclei has been attempted by a rela- 
tivistic field theory for the nuclear many-body problem. 
Reviews of this work have been given by Serot & Walecka 
(1986) and Gambhh et al. (1990). 

Relativistic Mean Field Theory starts from an effec- 
tive Lagrangian containing the nucleonic and mesonic de- 
grees of freedom. It is a relativistic analogue of the con- 
cepts of density dependent Hartree-Fock calculations with 
Skyrme forces. Because of its proper treatment of the spin- 
orbit splitting this model is expected to be more reliable 
than the non-relativistic models in predicting yet unknown 
properties of nuclei far from stability which are impor- 
tant in astrophysical situations. Also in some respects this 
method is simpler than the Skyrme type calculations, since 
the RMF method involves only local quantities such as lo- 
cal densities and fields. The binding energies and nuclear 
charge radii calculated by these methods agree well with 
the experimental values (see e.g. Table m. The density 
distributions of doubly magic spherical nuclei also agree 
well with electron scattering data. In addition, apart from 
the ground state properties of spherical nuclei, RMF re- 
produces the right ordering of the single particle spectra 
in adjacent odd mass nuclei, without any additional pa- 
rameters. It also describes in a quantitative way nuclear 
deformations and superdeformations. 

The microscopic RMF approach is used here to cal- 
culate nuclear properties at zero temperature, both for 
isolated nuclei, an approximation which holds well at the 
end of the Si-burning stage as well as for the case with 
higher density but still zero temperature where the nuclei 
can be assumed to be distributed in a lattice, and where 
modifications due to Lattice Coulomb energy have to be 
taken into account. 

Within relativistic mean field theory (Sheikh et al. 
1993) the nucleus is described a an ensemble of nucleons 
described by Dirac spinors and moving independently in 
self-consistently determined meson- and electro-magnetic 
fields. The Lagrangian, describing this system is given by 



L = Lb 
with 



C 



M 



c 



BMi 



Cb = ijji{i"f^dfj, -mN)^i 



(11) 



(12) 



-R^'^Rfj.i, 



\^l~P''~Pt.-\F^''F^.., 



— — 1 - T3 



(13) 



(14) 



mjq is the nucleon mass, m„, m^, and nip are the masses 
of the CT-, w- and p-mesons, g^ gui and gp are the cor- 
responding meson-nucleon coupling constants. Tilde indi- 
cate the isospin degree of freedom and fJ^", Rp^, Fp^ are 
the field tensors. U{a) is the non-linear potential for the 
(T mesons. It takes into account the density dependence in 
a phenomenological way and it has the form (Boguta & 
Bodmer 1977): 



U{a) ^ \m-y 



1 3 



1 4 

-^gzo 



(15) 



£ 



1 



M 



-d^^adpG - U{a) ~ -n^^^^p, + -Tv^'^oj^^u,^ 



For the three mass parameters mi^^ m^ and rUp entering 
the Lagrangian we use the experimentally known values. 
In general, the other parameters are obtained by fitting 
ground-state binding energies and charge radii of a few 
spherical nuclei. Several sets of such parameters are given 
in the literature. The parameter set NLl (Reinhard et al. 
1986) has turned out to be very successful in the valley of 
/3-stability. Here we use the parameter set NL-SH, which 
was obtained by fitting the binding energy, the charge- 
radii and also the neutron-skin radii of a few spherical 
nuclei (Sharma et al. 1993). In table Qwe show that both 
the NLl and the NL-SH parameter sets reproduce well 
the experimentally known nuclear charge-radii (Nadjakov 
et al. 1994) of Ni isotopes. We also compare (see Table |7|) 
with results obtained from the more recently determined 
parameter set NL3 (Lalazissis et al. 1997), which is often 
used in the literature for the study of exotic nuclei. 

While the set NL-SH reproduces well nuclear binding 
energies, we find that the incompressibility coefficient Kq 
of bulk nuclear matter as predicted by the set NL-SH is 
much larger than the value adopted earlier, where a value 
of Ko = 180 MeV has been used e.g. in BCK. 

A comparison of the neutron chemical potentials /i„ 
with the results of the compressible liquid drop model 
suggests that the RMF results are not very sensitive to 
the value of the adopted compressibility coefficient, (see 
Sect.|). 

However, it is to be noted that the results presented 
in this work are meant to apply to physical conditions 
prevalent in the collapsing core up to p = 10^^ g/cm'^, i.e. 
well below the density at which the core reaches maximum 
compression and rebounds. Since the nuclear composition 
of the core at this stage can still be treated as consisting 
of individual nuclei plus Coulomb lattice corrections, one 
may use a relatively high value of bulk compressibility as 
in NL-SH without worrying too much about the effects 
that it may have on the strength of the shock generated 
by the re-bouncing core. 
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Since mean-field calculations take into account only 
the long-range part of the nucleon-nucleon interaction the 
short-range pairing correlations have to be incorporated 
in addition. In the present work, the simplistic monopole 
pairing interaction is used and the resulting equations are 
solved in the BCS approximation. Because the value of 
the coupling constant in the pairing interaction is not 
known precisely, pair-gaps obtained from odd-even mass 
differences are generally used as a measure of the pairing 
strength. Initially, the calculations were started off using 
pair-gaps obtained by the 4-point method of Bohr & Mot- 
telson (1969) from the experimental values of the bind- 
ing energies, wherever known. However, the experimental 
binding energies are not known for most of the nuclei of 
interest and we have adopted the pair-gaps suggested by 
Moller et al. (1995) using isospin asymmetry arguments. 
For the cases where the experimental binding energies are 
known, these and the theoretically computed binding en- 
ergies are displayed in Fig.^j. 

In cases where multiple energy solutions were obtained, 
the solution corresponding to the lowest energy was taken 
as our preferred solution. It is found that deformations of 
the lowest energy solution predicted by the RMF calcula- 
tions are similar to those predicted by the Finite Range 
Droplet Model (FRDM) in Moller et al. (1995), which 
leads in this case practically to the same results as the 
Finite Range Liquid Drop Model (FRLDM). 

The values of the pairing gaps A„ and Ap used in 
the entire series of calculations have been recorded in ta- 
bles ^, ^ ^ and 13 for a series of isotopes ranging from 
^'*Mn25 to Ge32. We also show in these tables the cal- 
culated values of the quadrupole deformation beta /3, the 
binding energies BEhmf and the chemical potentials ^„ 
and Up. The binding energies are compared with the ex- 
perimental values BEexp. The data for the Ni isotopes is 
especially useful because these isotopes are spherical and 
in later sections we will compare the RMF calculations 
with macroscopic models as the macroscopic liquid drop 
model (BBAL, BBCW), which do not take into account 
effects due to deformation. 

In Table ^ we present for the chain of Ni- isotopes calcu- 
lations with the parameter set NL3 Lalazissis (1997). For 
most of the elements in this chain the binding energies 
obtained with the new parameter set NL3 are somewhat 
better, the set NL-SH however seems to be superior for 
very large neutron excess. Comparing with table ^ we find 
for the chemical potentials relatively small differences of 
around 100-200 keV between the two parameter sets. We 
therefore do not expect dramatic changes in using the pa- 
rameter set NL3 as compared to NL-SH, which is applied 
in all the other investigations in this work. 

How was this range of isotopes chosen? The significant 
difference of nuclei in the collapsing core from those in the 
laboratory is that the former are more neutron-rich. The 
typical electron fraction Ye ranges from 0.42 to 0.38 in 
the early stages of the collapse rather than from 0.46 to 



0.5 as is the case with the more symmetric laboratory nu- 
clei. Nevertheless they are still far less neutron-rich than 
matter in neutron stars where the Yg could be near 0.05. 
The nuclear statistical equilibrium gives rise to a fairly 
broad range of nuclei at any given instant of time, and 
in a given mass-zone of the collapsing core which can af- 
fect the (potentially) observable neutrino spectrum. The 
range of nuclei chosen here corresponds to values of the 
lepton fraction Yj, which for this pre-trapping range of 
temperatures and densities, equals the electron fraction 
Ye, ranging from 0.42 to 0.38. They are all mid-/ — p shell 
nuclei which have been shown (Hix & Thielemann 1996) 
to dominate the core nuclear configuration for this stage 
of the stellar evolution. 

4. Comparison of Chemical Potentials 

If a detailed simulation of the stellar collapse is carried out, 
where electron-capture and beta decay rates of individual 
nuclei have to be taken into account in a network calcula- 
tion, the EOS based on macroscopic spherical liquid-drop 
models cannot be expected to reproduce accurately (i.e. 
with a deviation of approximately 1 MeV) the nuclear 
chemical potentials (/Lt„, /Xp and (1) and binding energies 
(the Wm values) in the low-density region. In the follow- 
ing sections we investigate the amount of deviation that 
can be expected due to shell effects and pairing corre- 
lations on nuclear properties under stellar conditions. We 
also discuss the possible improvements which can be made 
on Wn using more advanced versions like the macroscopic- 
microscopic mass formulae, e.g. the Finite Range Liquid 
Droplet Model (FRLDM) of Moller et al. (1995), which 
account for deformation effects. 

The RMF computations presented in Sect. ^ are com- 
pared here with values of /j,„ and jl calculated from the 
models discussed in the previous sections, extended to the 
case of isolated nuclei, i.e. ignoring the lattice corrections. 
Comparison is also made between RMF calculations and 
the finite temperature BCK-EOS, at low temperatures 
and densities (pio = 0.1 and T = 0.1 MeV) using both 
the values of K^ and Sy adopted by BCK in their cal- 
culations and the more recent values used by Sharma et 
al. (1993). These comparisons show that the binding ener- 
gies as well as the neutron and proton chemical potentials 
have considerable differences from the values predicted by 
the macroscopic models of BBAL. The results from RMF 
calculations are in better agreement with the predictions 
from FRLDM calculations that are based on more recent 
compilations of the laboratory data on nuclei. 

4.1. Comparisons with the Spherical Liquid-Drop Models 

It has been stated earlier that the relations for /x„ and /i 
in Eqs. ^ and (^) are applicable for nearly symmetric, 
incompressible nuclear matter, confined by spherical sur- 
faces. The bulk energy term per nucleon Wbuik was taken 
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to be dependent on the nuclear volume energy and on the 
nuclear asymmetry energy (see BBCW Eqs 2.1 to 2.10, 
especially 2.4) given by Eq. |^. A derivation of fi„ using 
this form of the bulk energy gives the following relation 
for the neutron chemical potential of a mean nucleus of 
mass no. A :- 



H„ = -16 + 29.3(1 -2j;)2 + 117.2 a;(l- 2a;) 
7x — 3 

I ^* sizex-^ ) 777^ \" — f^nlvol i f^nlsize^ 

3(1 — X) 



(16) 



where Wsizeix) — 75.4a;^(l — x)^^^ is the combined sur- 
face and Coulomb energy (without lattice correction, i.e. 
fp[u] — 1) for the "mean nucleus" and po{x) is the nu- 
clear density. Fig^ displays the chemical potential ob- 
tained form Eq. (JK) for Ni isotopes, alongside with their 
RMF and BCK-EOS counterparts. It should be noted here 
that because of an error in the surface terms, the /z„ from 
the BBAL equations do not agree well with the RMF re- 
sults (in fact, they are off by up to 2 MeV). 

The RMF results are also compared with the BCK- 
EOS in Fig.g. For generating this set of numbers for the 
BCK-EOS, the density and the temperature were set at 
10^ g/cm^ and and 0.1 MeV so that we are as close as pos- 
sible to the end-stages of Si-burning, with their nuclear 
asymmetry energy coefhcient taken as Es = 31.5 MeV 
and the incompressibility as K = 180. However, the RMF 
calculations correspond to a incompressibility coefficient 
K = 354.95 and Sy = 36.1 for reasons discussed in Sect. 
H. Values for the nuclear incompressibility deduced from 
experiments have shown a large variation, Kq ranges from 
180 to 300 MeV. Since the objective of the comparisons 
here is to evaluate the deviation that can be expected in 
the LDM-EOS (compressible and incompressible) due to 
exclusion of shell and pairing effects, the //„ and jl from 
BCK-EOS were recalculated using different values of Sy 
and Kq. It was found that the results were more sensitive 
to a variation in Sy than to Kq. In order to obtain a suit- 
able value of Sy (used in BCK-EOS) which would better 
reproduce the RMF chemical potentials, a fit was made of 
the following expression of /IIrmf- 



fl = ASy{l - 2x) - 290 



x^{l-x)^ 



2 l-2x 



X 1 



(17) 



This gives an "average" value of Sy (over the range of 
isotopes for Mn, Fe, Co, Ni, Cu, Zn, Ga, Ge considered in 
the tables |, |, |, ^ of 30.4 MeV. It is found that using 
the low value oi K = 180 MeV, the "best fitted" value of 
Sy is 30.34 MeV. Figs.|| and ^ compare ^„ and p. for Ni 
isotopes using this parameter set in BCK-EOS with RMF 
results. 

4.2. The Finite Range Liquid Drop Model (FRLDM) 

The various forms for the total binding energy of a nu- 
cleus Wn that have been quoted in Sect. have all been 



for spherical nuclei without shell and pairing corrections. 
A refinement of these simple liquid drop mass formulae 
of Myers & Swiatecki (1966) and Myers et al. (1995) has 
been made using the microscopic-macroscopic methods. In 
this approach, the total nuclear potential energy is taken 
to be the sum of the energy of a term which deals with the 
bulk properties of nuclear matter and a set of terms which 
deal with both the shell and pairing effects as well as ef- 
fects due to nuclear deformation. Here the macroscopic 
Finite Range Liquid Drop Model FRLDM is used to cal- 
culate neutron chemical potentials /z„ in the stellar core 
for the same range of nuclei for which RMF calculations 
were made in the previous section. This approach is jus- 
tifiable for pre-trapping densities where nuclear matter is 
treated as an ensemble of isolated cold nuclei with a small 
drip neutron fraction. At higher densities two additional 
effects would have to be taken into account: (a) the drip 
neutron contribution as given in Eq. |lj and (b) the lattice 
corrections made in the Coulomb term. The Finite Range 
Liquid Drop model is an extension of the LDM of Myers 
& Swiatecki (1966) with modifications made for the finite 
range of the nuclear forces by means of a folded Yukawa 
potential with exponential term. 

Recently there has been a further modification of the 
FRLDM (Myers & Swiatecki 1995) to the FRDM, the Fi- 
nite Range Droplet Model (MoUer et al. 1995). This has 
been done in order to account, among other things, for the 
corrections due to Coulomb redistribution effects and the 
effect of nuclear incompressibility on neutron and proton 
radii. This improves the predictions of higher moments of 
nuclear deformation and binding energies in the FRDM 
model. However, it also depends on a larger number of 
parameters. Hence, the older FRLDM model is used for 
in most of our calculations. From Fig. M and Fig. @ we see 
that these models make practically no differences for the 
present calculations. 

In the FRLDM model, the macroscopic component of 
the total nuclear binding energy Wn = WmaciZ,N) is 
given by 

WrnaciZ, N) = -a„(l - Kyl^) + as(l - At,/2)BiA2/3 

XjA Z=N, odd 
otherwise 



-Ca(^-Z)-f M/ |/| + 



where the average pairing energy Apair is obtained from 
nuclear binding energies using the 4-point method of Bohr 
et al. (Sect. |h and is given by: 



(18) 



where / = {N — Z)/A is the relative neutron- proton ex- 
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Wmac{Z,N) is the volume energy, the second the sur- 
face energy, with a correction for surface asymmetry, the 
third term is due to the Coulomb energy, the fifth is the 
Coulomb exchange term. The other terms are (in that or- 
der) due to proton form factor correction to the Coulomb 
energy, the charge asymmetry energy, and the Wigner 
term respectively. The numerical values of the various pa- 
rameters used in this model are quoted in table ||. The 
proton form factor correction to the Coulomb energy is 
given by 
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(19) 



Both coefficients Bi and B^ are shape dependent quanti- 
ties. The factor Bi is the relative generalized surface or 
nuclear energy. It takes into account the finite range a of 
the nuclear forces and is given by: 



Bi = 



87rV2a4 
2- ■ 



r-ri \ exp(-|r-r'|/aj3^,3 / 



r'l/a 



-d'rd'r' 



(20) 
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l-(l + i'^-^'' 



2 ade 



and B^ is the relative Coulomb energy and is given by: 

-1 = 4-5/3 

B,= ^4^-^x (21) 

j cxp(-|r-r'|/a<ien) 

Using Eq. (|^), we find the following expression for the 
volume contribution to /i„ : 

fJ-nlvoi = -16 + 123.04(0.5 - xf + 246.08a;(0.5 - x) (22) 

Note that the bilinear term in this expression reduces ap- 
proximately to the linear term in Eq. (R) when symmetric 
nuclei are considered {x ~ 0.5). The surface etc. contri- 
butions to Unltotai were obtained by using Eq. (||) on the 
surface energy, the Coulomb energy and the Coulomb ex- 
change correction components of the FRLDM expression 
for total energy quoted above. The small corrections due 
to the shape dependent terms Bi and B3 were ignored, as 
were the other terms in the expression for the total nuclear 
binding energy. 

In Fig.0 the functional form of fi„ in Eq- (P) using 
FRLDM is compared with the corresponding RMF values 
for the spherical Ni isotopes and in Figs.0, M and nu for 
the deformed Mn, Fe, and Ga isotopes. It is found that the 
FRLDM model reproduces the RMF values better than 
either the BCK-EOS (with the "best" parameters em- 
ployed by them) or the spherical-LDM results in Eq. (|^) . A 
comparison of the binding energy predictions of both the 
spherical, incompressible LDM-EOS and FRLDM model 



is shown in Figs.O and n2[ The latter models reproduce 
the binding energies with an accuracy of ~ 1 MeV. 

In Fig.|| we compare the chemical potential of neu- 
trons in the Ni-isotopes calculated within RMF-theory us- 
ing parameter NL-SH with the results of the Finite Range 
Droplet Model (FRDM) of MoUer et al. (1995). 



5. Equation of State at Finite Temperatures 

At low densities, the neutrons, protons and nuclei can 
be described by an ideal Boltzmann gas. The supernova 
matter is then in nuclear statistical equilibrium at fixed 
electron concentration Yg, where the electrons, positrons 
and neutrinos (when they are trapped) are non-interacting 
Fermions (El Eid & Hillebrandt 1980). This approach is 
justified as long as the nucleon vapor is dilute, i.e. up to 
about a tenth of the nuclear saturation density. It does 
not strongly modify the nuclear properties. Higher den- 
sities and temperatures require an alternate microscopic 
description which accounts for interactions between all the 
nucleons present. Approaches such as equilibrium in the 
bulk, the compressible liquid drop model, the Thomas- 
Fermi, or the Hartree-Fock method have been used in the 
literature. 

In the macroscopic approach, the density of a given 
species is obtained from equations of the Saha-type where 
one requires partition functions which properly account 
for the internal excitations of the nuclear species. One also 
needs the nuclear ground state energies, i.e. the binding 
energies of these nuclei. Usually one uses semi-empirical 
mass formulae for the binding energies and Fermi gas 
based level densities for the nuclear partition functions. 
The quality of these ingredients determine the level of ac- 
curacy of the EOS and its range of applicability. While the 
macroscopic method may be adequate for the computation 
of stellar collapse dynamics, one has to take into account 
shell-effects and pairing correlations in order to calculate 
threshold Q- values for electron captures on nuclei, if one is 
interested in the determination of the neutrino spectrum. 
In addition, ideally one would like to have an approach 
which is applicable in both low and high density situations 
encountered in supernova physics. From this point of view 
a microscopic approach like the RMF method, which sat- 
isfies both criteria, is a preferred method. It can also be 
used to compare results obtained from other microscopic 
methods. 

It can be shown that in the case of low entropy of the 
stellar core matter, e.g. at the stage of collapse, a multi- 
phase equilibrium such as between nuclei and nucleonic 
vapor can be conveniently studied by using a low tem- 
perature expansion for the dense nuclei, along with a high 
temperature expansion for the nucleonic vapor (Vautherin 
1994). Vautherin's description of the problem is a mean- 
field approximation and a simple extension of our RMF 
calculations to finite temperatures is thus possible. 
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For a spin saturated, symmetric nucleus (N=Z) the 
mean field potential from a Skyrme force with a zero range 
two-body attraction and a three-body repulsion, with co- 
efHcients to = -983.4 MeV fm^ and ig = 13105.8 MeV 
fm^ is given by (Ring & Schuck 1980): 



C/(r) = -topir) 



16 



hp\r). 



(23) 



With this potential the single particle orbitals 0j(i") ^^'^ 
given by the solution of the Schrodinger equation with 
eigenvalues e^. The Fermi occupation numbers fi at finite 
temperatures are given by /j = (1 + cxp((ei — p)/kT))^^ 
and the chemical potential /i is fixed from the particle 
number ^^ /j = A. This gives the self consistent nucleon 
density at a finite temperature as p(r) = Ylifi \'t'i{^)\ ■ 
The entropy of the nucleus obtained in a similar way by 
the expression for independent fermions: 



S = - A: ^{/, log /, + (l-/,)log(l -/,)}. 



(24) 



For the uniform nuclear matter approximation, the density 
of the nucleons and kinetic energy densities are found as 



and 

r = 47rgA-^F5/2(/?C/-a) 



(25) 



(26) 



where g=4 is the spin and isospin degeneracy factor 
(2S'+1)(2/+1), A = \2TTri^/mkTY/^) is the thermal wave 
length, [3 = l/kT, and a = (3p. The Fermi-integral Fn is 
given by: 

/>oo 

F„(x) = (2/0F) / y"-i(l + exp(y + x))-^dy. (27) 

Jo 

Once the density and the temperature of stellar matter 
is specified, the argument {j3U — a) of the Fermi func- 
tion is determined, which in turn can be used to compute 
kinetic energy density, entropy density, pressure, and the 
free energy density of the system. 

A consideration of the phase diagram, i.e. the isother- 
mal lines in P — p plane, indicates that at temperatures 
below a critical temperature a configuration with constant 
density does not correspond to the lowest value of the free 
energy for some density ranges. One has a two phase sys- 
tem consisting of a uniform density nucleus and the nu- 
cleonic vapor in equilibrium with each other as a more 
favorable state. Under such conditions, one has: 

p{pN,T)^p{p,,T) and P{pN,T) ^ P{p,,T). (28) 

So far, two effects which are important for nuclei have 
been left out, the finite size (or surface) effects and the 
Coulomb interactions. Levit and collaborators (Levit & 
Bonche 1985; Besprosvany & Levit 1989) have considered 
the nucleus as a homogeneous system oi A = N + Z nu- 
cleons inside a spherical volume fl of radius R which is 



in equilibrium with the nucleonic vapor of volume Qy and 
density p^. Because of the large incompressibility of nu- 
clear matter it can be shown that the change in the den- 
sity of the nuclear phase caused by Coulomb and surface 
tension effects is rather small. Therefore in the equilib- 
rium equations between the nuclear and vapor phases, it 
is adequate to evaluate the Coulomb and surface terms 
for the pressure and the chemical potentials at saturation 
density po and R = Rq = (S/Anpo)^/^ (Vautherin 1994). 
Using the Gibbs-Duhem relation pdp/dp = dP/dp one 
can also eliminate the density of the nucleus phase in the 
vapor-nucleus equilibrium equations, and observing that 
the coefficients oi p ~ po in the expansion of the nuclear 
pressure and of the chemical potential have to be equal. 
Thus one has: 



Eo ^^{kTf 



1 f2a{T) 
Ro V Po 



^4~ 



(29) 



The RMF calculations outlined in the previous sections 
automatically account for the surface and Coulomb effects 
in the nuclear chemical potential calculations (such as the 
last term within parentheses in the above equation) . Since 
Tc is quite high, it is adequate to use in the computation of 
the equation of state at low densities the zero temperature 
surface tension effects already built into the RMF code 
for the nuclei, i.e. a{T) = a(0). Therefore the only finite 
temperature correction is the second term on the RHS of 
the above equation. 

5.1. Lattice Size Effects 

As an attempt to extend the RMF results presented in 
section ||, to higher densities, the density dependent lat- 
tice correction should be incorporated into the Laplace 
equation for the electric potential: 



-AAO(r)=pp(r)-pe(r) 



(30) 



where A" is the Coulomb potential, Pp(r) is the proton 
charge density and yOe(r) is the electron charge distribu- 
tion density in the Wigner-Seitz cell. Here, the lattice cor- 
rection has been incorporated into the Coulomb energy 
part of the total energy by the inclusion of the multiplica- 
tive factor Fiattice = (1 - |(^)^/^ + 5(^)) where p and 
Ps are the stellar matter and nuclear saturation densities 
respectively. In Fig.O we show the values of p for Ni iso- 
topes along with the results from Eqs. (Q) and (||) for the 
matter density p = 10"'^^ gcm"'^. 

At densities near « 10^^ g cm^'^ the distance between 
individual nuclei becomes so short that the nuclear mo- 
tions begin to be correlated due to the Coulomb interac- 
tion between the nuclei and at somewhat higher densities 
the nuclei arrange themselves on the lattice sites of a crys- 
tal. Thus, towards the end of the density regime that one 
would be interested in for computing early (un-scattered) 
neutrino spectra, EOS calculations have to take into ac- 
count screened Coulomb interactions. This is usually done 
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by calculating thermodynamic quantities in the Wigner- 
Seitz approximation. Stellar matter with a plasma param- 
eter r = Z^e^/RkT > 155, arranges itself in a lattice and 
the Coulomb interaction energy can be written at zero 
temperature as: 

3"! 



Er^ 



3ZV 
5 R 
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3R 
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(31) 



This therefore motivates the inclusion of the density de- 
pendent Fiattice correction in the RMF equations. 

5.2. Asymmetric Nuclear Matter and Approximate 
Solution of the Equilibrium Equations 

To describe nuclei with different neutron and proton num- 
bers as in the neutron-rich situation encountered in the 
stellar collapse one has to include the bulk asymmetry 
term arxj^ in the mass formula, where xn = [N — Z)/A is 
the nuclear asymmetry parameter. By writing the equilib- 
rium equations explicitly in terms of low and high temper- 
ature expansions in the nuclear and in the vapor nucleonic 
chemical potentials, one is able to relate the difference in 
the nuclear component of the neutron and proton chemical 
potentials to the asymmetry parameter xn as: 

y.Nn- tJ'Np^'^ar{T)xN, (32) 

where, ar(T) is the second derivative of the free energy 
with respect to xn'. 

, , 1 I , , , 1 , vr^ (fcT)2 ^ ^ 

ar{T) = -TN--t^{x^ + l/2)pN-—t^pl + — ^-j^{ii) 

where xq — 0.48. Similarly the vapor asymmetry x^ and 
nuclear asymmetry can be related as: 

Xt, — tanh(2aT-a:^Ar//cT), (34) 

while the vapor density can be related to the average nu- 
clear chemical potential /zat = [pNn + f-Np)/2: 

(35) 



p, = g\~^{l - xl)-^/^expifiM/kT). 

Here fj.]\[ defined above can be approximated to the first 
order in xn as: 
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A'-:^^-^U^-^.(36) 
Po V 2Tf J 12Tf 

It has been shown (Vautherin 1994) that it is a good ap- 
proximation to replace Or by its value at zero temperature 
and normal density, since the dependence of a^ on tem- 
perature is weak. The pressure in the vapor phase is given 
by the high temperature approximation: 



Pv ^ PvTil 



45V2 



(!+<) + 



+ iioP.'(3 - xl{2xo + 1)) + \t^pl{l - xl). (37) 

Although the total pressure is dominated by the lepton 
pressure, one may estimate the nucleonic component from 
the above expression. 



6. Summary 

Using the relativistic mean field method, we have evalu- 
ated in this paper properties of individual nuclear species 
which play a significant role for the physical and ther- 
modynamic properties as well as the dynamics of a pre- 
supernova collapsing stellar core. In particular, we have 
determined the nuclear chemical potentials p„ and pp, 
which, along with the nuclear weak interaction strengths 
control the neutronization of the matter. We have focussed 
our interest on properties of nuclei present in the core at 
a density range 10^ to 10^^ g/cm^ which corresponds to 
the temperature range of ~ 0.2 to 1.5 MeV. In this range 
the core composition is dominated by heavy, neutron-rich, 
f — p shell nuclei in a Saha equilibrium, with a low density 
vapor of free neutrons, and an even lower density vapor 
of free protons. In this temperature and density range, 
matter can be treated as an ensemble of non-interacting 
Boltzmann particles, at least at pio — 0.1 — 100 gjcrrt" . 
This implies that nuclear properties under these condi- 
tions are determined by those of individual nuclei, and 
structure effects like shell effects, pairing and deformation 
play a significant role. 

Our study has been motivated by the fact that the 
existing equations of state used for numerical studies of 
stellar core collapse are based on macroscopic liquid-drop 
models, which generally do not take into account shell, 
pairing and other nuclear structure effects - which play an 
important role at lower energies. These EOS arc suitable 
for matter at T > 1.5 MeV, and p > 10^^ g'/cm^, when 
the nuclear matter undergoes a phase transition from 
the state where individual nuclei are present in a gas of 
dripped nucleons, to the state where the nuclei merge into 
infinite nuclear matter (see Sect. and (BBAL, BBCW, 
Cooperstein (1985)). Because shell and pairing effects are 
already washed out at temperatures beyond 1 MeV and at 
densities beyond trapping density i.e. 10^^ g/cm'^, ideally, 
a low-density, zero temperature EOS would be desirable, 
which smoothly goes over to the high density, high tem- 
perature equations of state based on the liquid drop model 
(see Swesty et al. 1994). 

In order to account for the shell and pairing effects per- 
sistent at low temperatures we studied the nuclear struc- 
ture effects in the framework of relativistic mean field the- 
ory, using a Lagrangian with non-linear self-interactions 
for the cr-meson. This method has been shown to provide 
a good description of the ground state properties of nuclei 
(Gambhir et al. 1990). The nuclei of interest range from 
s^Mnss to **iGe32 (see Tables |, |, |, and |), predomi- 
nantly neutron-rich f — p shell nuclei. It is the first time 
that the RMF method is used to calculate the binding 
energies of nuclei with a neutron excess in the range of 
interest. 

In order to test the validity of high-density equations 
of state at low densities BBAL we compare the chemical 
potentials pn and p with those obtained in RMF theory 
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and we find considerable deviations. To study the influ- 
ence of nuclear compressibility we have also compared the 
chemical potentials HuIrmf and /ip|_RMF with those of 
the compressible liquid-drop model based on the EOS of 
BCK. To test the effects of a variation in the nuclear in- 
compressibility Kq and in the volume symmetry energy 
Sv , we computed the results of BCK with a range of values 
of Kq and Sy. We find that the results are more sensitive 
to the variation in Sy than in Kq. From a fit of jj-Irmf to 
the analytical expression for fi\BCKi we find an "averaged" 
value of Sv of 30.34. Further, to account for nuclear defor- 
mation effects we use the FRLDM expression for the nu- 
clear binding energy to compute /i„ and /t. A comparison 
of RMF values with this "extended" macroscopic model 
demonstrates the influence of shell and pairing effects on 
fin and on fi (see Figsjll| and |lj). 

Finally, in order to take into account the first order 
Coulomb correlations between the nuclei at high densities, 
we include a lattice correction in our RMF computations 
of fi. A comparison of these lattice corrected RMF results 
with the BBAL/FuUer EOS for fi is shown in Figjl3[ We 
have also outlined the reasons for obtaining the relatively 
small corrections that are needed for finite temperature 
corrections to the EOS at low entropies relevant for the 
early stages of a massive star's core collapse. Methods to 
generalize the results to finite temperature EOS calcula- 
tions, using the chemical potentials and other quantities 
obtained from RMF calculations have been indicated. 

Many groups have computed and analysed the collapse 
of the stellar core of a massive star to evaluate the effects of 
various physical processes in determining conditions prior 
to core-bounce. As collapse is shown to be a nearly ho- 
mologous one (see BBAL and references therein), a one- 
zone collapse is adequate for the investigation of ther- 
modynamic and nuclear properties (see e.g. EH (1980); 
HiUebrandt et al. (1984) which used the EOS of Wolff & 
Hillebrandt (1983); Ray et al. (1984); Fuller (1982); etc). 
The EOS of dense stellar matter affects the electron cap- 
ture rate and the evolution of the lepton fraction mainly 
through the influence of the chemical potentials of the 
free nucleons (see Murphy (1980); Sutaria & Ray (1997)). 
For example, the thermodynamic quantities at the begin- 
ning of stellar collapse for the RMF results and those of 
EH are compared in Table H. However, the intrinsic elec- 
tron capture rates on nuclei employed by different authors 
have been different and have evolved considerably over the 
years. It is therefore difficult to extract the purely EOS 
related effects on neutronisation of different calculations. 
The stellar collapse calculation using the full RMF EOS in 
a self-consistent manner with the latest electron capture 
rates is in preparation and the details of the evolution of 
nuclear and thermodynamic variables will be published 
elsewhere. 
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Table 1. The nuclear charge radii for Ni isotopes. The 
table gives the experimental values and those computed 
in this work from RMF calculations with the parameter 
sets NLl and NL-SH. 



F. K. Sutaria et al.: Nuclear properties in stellar collapse: RMF 

The parameter set for macroscopic FRLDM 



Table 

model 



Quantity Brief definition 



Value (Unit) 



Neutron 


Experimental 


RMF 


RMF 


number 


Nadjakov et al. (1994) 


(NLl) 


(NL-SH) 


30 


3.7827 ± 0.0036 


3.804 


3.791 


32 


3.8177 ± 0.0047 


3.797 


3.794 


33 


3.8283 ± 0.0030 


3.805 


3.802 


34 


3.8475 ± 0.0050 


3.812 


3.809 


36 


3.8678 ± 0.0051 


3.830 


3.826 



Reinhard P.-G., Rufa M., Maruhn J., Greiner W., Friedrich J., 

1986, Z. Phys. A323, 13. 
Ring, P., Schuck, P., The Nuclear Many-body Problem, 

Springer Verlag New York, 1980. 
Serot B. D., Walecka J. D., 1986, Adv. Nucl. Phys. 16, 1. 
Sharma M.. M., Nagrajan M. A., Ring P., 1993, Phys. Lett. 

B312, 377 
Sheikh J. A., Maharana J.P., Gambhir Y.K., 1993, Phys. Rev. 
C48, 192. 

, Ph.D. thesis, Mumbai University 1997a. 
, Ray A., 1997, Phys. Rev. Lett. 79, 1599 
Sheikh J. A., Ray A., 1997b, Nucl. Phys. A621, 



Lattimer, J.M., Myra, E.S., 1994, Ap. J. 425, 



Sutaria F. K., 
Sutaria F. K., 
Sutaria F. K., 

375c. 
Swesty, F.D., 

195. 
von Groote, H., Hilf, E. R., Takahashi, K., 1976, At. Data Nucl. 

Data Tables 17, 418. 
Vautherin D., 1994, in Supernova, Les Houches Session LIV, 

ed. S. Bludman, R. Mochkovitch and J. Zinn- Justin 
Wolff, R. G., Hillebrandt, W., 1983 in "Nucleosynthesis", eds 

Arnett and Truran, Univ of Chicago Press, p. 131. 



as 

Ks 

ao 

Ca 

ro 
rp 
a 



Ode 



Volume energy constant. 

Volume asymmetry constant. 

Surface energy constant. 

Surface asymmetry constant. 

A*^ constant. 

Charge asymmetry constant. 

Nuclear-radius constant. 

Proton r.m.s. radius. 

Range of 

Yukawa- plus- exp onent ial 

potential. 

Range of Yukawa 

function for generating 

charge distribution. 

square of the e~-charge 



16.00126 MeV 

1.9224 MeV 

21.18466 MeV 

2.345 MeV 

2.615 MeV 

0.10289 MeV 

1.16 fm 

0.80 fm 

0.68 fm 



0.70 fm 



1.4399764 MeV fm 
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Table 3. Results of RMF calculations for various isotopes 
of Mn and Fe, giving the values of /x„, /ip and the defor- 
mation parameter 13. Also tabulated are the values of the 
fixed neutron and proton pairing-gaps (A„ and Ap) used 
in these calculations. The last two columns display the ex- 
perimental and calculated values of the binding energies 
{BE^xp. and BE^mf respectively) 



Isotope 


A„ 


A^ 


Hn 


Mp 


13 


BEexp. 


BEkmf 


"''Mnas 


1.00 


1.50 


-8.6790 


-8.6260 


0.1265 


-471.84940 


-471.64801 


'^''Mnas 


1.00 


1.50 


-8.7800 


-9.2290 


0.2082 


-482.07599 


-480.91101 


^''Mnas 


1.20 


1.50 


-7.9420 


-9.9240 


0.2187 


-489.34659 


-488.92999 


''^Mnas 


2.00 


1.50 


-7.9260 


-10.5230 


0.2091 


-497.99701 


-497.51801 


""Mnas 


1.20 


1.50 


-7.3470 


-11.4300 


0.2213 


-504.41000 


-503.65399 


'''Mnas 


2.00 


1.25 


-7.1270 


-11.9500 


0.2130 


-512.12903 


-511.54700 


'°Mn25 


1.20 


1.25 


-6.2840 


-12.8120 


0.2110 


-517.62000 


-516.78497 


'^'Mnas 


2.00 


1.25 


-6.3820 


-13.1600 


0.1896 


- 


-524.057 


''"Mnas 


1.25 


1.25 


-6.0970 


-13.7840 


0.1354 


- 


-528.286 


''^Mnas 


1.50 


1.00 


-6.0110 


-14.3490 


0.1097 


- 


-534.347 


'^Mn25 


1.00 


0.50 


-5.8150 


-15.0530 


0.0830 


- 


-539.142 


''Mn25 


0.75 


0.50 


-5.2780 


-15.5570 


0.0458 


- 


-544.673 


""Fe26 


2.50 


1.50 


-12.0530 


-7.1700 


0.0003 


-471.76392 


-470.85400 


^^Fe26 


1.50 


2.50 


-10.5820 


-7.8950 


0.0000 


-481.06192 


-481.15302 


^•^Feae 


2.50 


2.50 


-9.5510 


-8.2400 


0.0000 


-492.25992 


-492.17700 


'^Fe26 


2.00 


2.50 


-8.9800 


-8.7980 


0.1531 


-499.90610 


-500.02802 


^«Fe26 


2.50 


2.50 


-8.7690 


-9.3750 


0.1514 


-509.95081 


-509.32599 


^^Fe26 


2.00 


2.00 


-8.1790 


-10.0260 


0.1915 


-516.53180 


-515.96600 


^°Fe26 


2.25 


2.00 


-7.8740 


-10.6300 


0.1839 


-525.34802 


-524.05402 


«^Fe26 


1.75 


2.00 


-7.2590 


-11.3220 


0.1853 


-530.93201 


-530.29102 


^^Fe26 


2.00 


2.00 


-7.1360 


-11.9780 


0.1421 


-538.97998 


-537.29602 


«^Fe26 


1.75 


1.75 


-7.0720 


-13.0450 


-0.0938 


-543.34998 


-542.30103 



14 



F. K. Sutaria et al.: Nuclear properties in stellar collapse: RMF 



Table 4. Results of RMF calculations for various iso- 
topes of Co and Ni, giving the values of /in, Mp ^n.d the 
deformation parameter (3. Also tabulated are the values of 
the fixed neutron and proton pairing-gaps (A„ and Ap) 
used in these calculations. The last two columns show the 
calculated and experimental values of binding energies. 



Isotope 


A„ 


A^ 


Hn 


Mp 





BEnMF 


B£'exp. 


''Co2r 


1.00 


1.50 


-9.2070 


-7.1380 


0.0003 


-486.91150 


-487.96399 


''Co2r 


1.02 


1.50 


-9.6860 


-7.2640 


0.1408 


-498.28741 


-497.01001 


^«Co2r 


1.17 


1.50 


-9.2300 


-7.7340 


0.1675 


-506.86050 


-506.35599 


5«Co2r 


2.00 


1.50 


-9.1760 


-8.4530 


0.1430 


-517.31409 


-516.29102 


''"Coar 


1.20 


1.25 


-8.6380 


-8.8190 


0.1820 


-524.80615 


-523.69702 


^^Coar 


2.00 


1.25 


-8.4150 


-9.6480 


0.1520 


-534.. 12738 


-532.99597 


«^Co2r 


1.15 


1.25 


-7.6400 


-10.1580 


0.1640 


-540.72498 


-539.52698 


«3Co2r 


2.00 


1.25 


-7.9050 


-11.4190 


0.0915 


-549.21198 


-548.28497 


«^Co2r 


1.20 


1.25 


-7.8090 


-12.4110 


0.0455 


-555.23602 


-554.36798 


«^Co2r 


1.50 


1.00 


-7.4920 


-13.2190 


0.0011 


-562.67999 


-562.22900 


^«Ni28 


2.50 


2.50 


-10.8710 


-5.7110 


0.0000 


-506.45901 


-507.02899 


^^Ni28 


2.00 


2.00 


-9.9250 


-6.3600 


0.0000 


-515.45892 


-514.91101 


«°Ni28 


2.50 


2.00 


-9.8810 


-6.9550 


0.0000 


-526.84735 


-525.61603 


"Ni28 


2.00 


2.00 


-9.3950 


-7.6600 


0.0000 


-534.66760 


-533.20502 


•''Ni28 


2.25 


1.75 


-9.2210 


-8.2490 


0.0000 


-545.26501 


-542.88800 


''^Ni28 


1.75 


1.50 


-8.9080 


-8.9340 


0.0000 


-552.10352 


-550.31299 


«*Ni28 


2.00 


1.50 


-8.6650 


-9.4930 


0.0000 


-561.76013 


-559.63000 


'''Ni28 


1.75 


1.50 


-8.3620 


-10.1220 


0.0000 


-567.85834 


-567.52698 


'''Ni28 


1.00 


1.00 


-8.0750 


-10.7160 


0.0000 


-576.83301 


-574.58398 


''Ni28 


1.50 


1.00 


-7.5840 


-11.2130 


0.0000 


-582.61902 


-583.01501 


IN 128 


0.75 


1.00 


-6.9360 


-11.7330 


0.0000 


-590.42999 


-589.70001 
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Table 5. Results of RMF calculations for various iso- 
topes of Cu and Zn, giving the values of /i„, fip and the 
deformation parameter (3. Also tabulated are the values of 
the fixed neutron and proton pairing-gaps (A„ and Ap) 
used in these calculations, as well as the experimental and 
calculated values of the binding energies. 



Isotope 


A„ 


A^ 


Hn 


Pp 


13 


BEexp. 


BEkmf 


'■■'CU29 


1.00 


1.50 


-14.0190 


-1.3950 


-0.0009 


-484.76001 


-485.70700 


^«CU29 


2.00 


1.75 


-11.0700 


-2.5080 


0.0037 


-497.11401 


-497.35001 


^«CU29 


2.00 


1.75 


-11.2610 


-4.0470 


0.1461 


-509.87601 


-508.19299 


«0CU29 


2.00 


1.75 


-10.6050 


-4.9190 


0.1650 


-519.93811 


-518.96899 


"CU29 


2.25 


1.75 


-10.5020 


-5.3060 


0.1394 


-531.64679 


-530.08899 


«^CU29 


1.75 


1.75 


-9.8260 


-6.3300 


0.1750 


-540.53400 


-538.87299 


«^CU29 


2.50 


1.50 


-9.6850 


-6.4610 


0.1374 


-551.38715 


-549.16901 


«*CU29 


2.00 


1.50 


-9.0460 


-7.2040 


0.1424 


-559.30316 


-556.85797 


«^CU29 


2.25 


1.50 


-9.1590 


-7.3710 


-0.1001 


-569.21216 


-567.38501 


«^CU29 


1.50 


1.50 


-8.8200 


-8.0180 


-0.0922 


-576.27814 


-574.91803 


«^CU29 


2.25 


1.25 


-8.5780 


-7.9580 


-0.0083 


-585.39697 


-583.46503 


""Znao 


2.50 


2.50 


-10.6040 


-5.8830 


0.1700 


-559.09912 


-557.58899 


'"Znso 


2.00 


2.50 


-9.8550 


-6.7130 


-0.1663 


-567.07898 


-566.32001 


"'Znso 


2.50 


2.50 


-9.7730 


-7.0710 


-0.1314 


-578.13885 


-576.87000 


67r7„ 

Znso 


2.00 


2.50 


-9.4760 


-7.5290 


-0.1161 


-585.19086 


-584.99200 


««Zn3o 


2.25 


2.25 


-9.2010 


-7.8250 


-0.0793 


-595.38898 


-593.99298 


^^Znso 


1.75 


2.00 


-8.7400 


-8.2670 


0.0159 


-601.87122 


-601.14203 


™Zn3o 


2.25 


2.00 


-8.4010 


-8.8290 


0.0169 


-611.08600 


-610.53101 


■"^Znao 


1.75 


2.00 


-7.8710 


-9.4300 


0.0090 


-616.91901 


-617.31897 


""Znso 


2.00 


1.50 


-7.6870 


-9.9300 


0.0764 


-625.80200 


-624.49799 


Znao 


1.25 


1.50 


-7.6010 


-11.0730 


0.0000 


-631.15002 


-630.94098 


J^Znao 


1.50 


1.50 


-7.2610 


-11.7630 


0.0000 


-639.51898 


-638.95001 


Znso 


1.25 


1.50 


-6.9250 


-12.4690 


0.0000 


-644.58002 


-645.87000 


"'Znso 


1.00 


1.00 


-6.4360 


-13.0960 


0.0001 


-652.41998 


-652.28497 


^^Znso 


0.20 


0.50 


-5.8320 


-13.7630 


0.1720 


-656.94000 


-657.80603 


'"Znso 


0.50 


0.50 


-5.4660 


-14.2850 


0.1497 


-664.06000 


-663.66101 
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Table 6. Results of RMF calculations for various iso- 
topes of Ga and Ge, giving the values of fin, fip and the 
deformation parameter (3. Also tabulated are the values of 
the fixed neutron and proton pairing-gaps (A„ and Ap) 
used in these calculations as well as the calculated and the 
experimental values of the binding energies, BE^xp. and 
^^'rmf respectively. 



Isotope Ari 



A„ 



pin 



fip 



BEe^ 



BEhmf 



""Gaai 


2.50 


2.50 


-9.7100 


-7.2800 


0.0897 


-601.99402 


-602.03601 


™Ga3i 


1.50 


2.50 


-9.2860 


-7.7540 


0.0010 


-609.64899 


-609.17798 


"'Gaai 


2.50 


2.50 


-9.0280 


-8.3030 


0.0148 


-618.95618 


-620.22803 


^^Ga3i 


1.50 


2.50 


-8.3770 


-8.8990 


0.0014 


-625.47729 


-626.45203 


^^Gaai 


2.50 


1.50 


-8.4390 


-9.4400 


0.1686 


-634.66302 


-635.96198 


^^Gaai 


1.50 


1.50 


-8.1890 


-10.0830 


0.1770 


-641.08002 


-641.78699 


"^Gaai 


2.00 


1.30 


-7.8210 


-10.5600 


0.1730 


-649.56598 


-650.75098 


™Ga3i 


1.00 


1.10 


-7.6520 


-11.0730 


0.1800 


-655.62000 


-656.83698 


''Ga3i 


0.75 


0.75 


-7.1730 


-11.6270 


0.1790 


-663.65997 


-664.02197 


'*Ga3i 


0.50 


0.50 


-6.7420 


-12.4590 


0.1680 


-668.88000 


-670.63800 


™Ga3i 


0.25 


0.25 


-6.0880 


-13.3830 


0.1540 


-676.14001 


-677.14899 


80 no 
Cja3i 


0.20 


0.20 


-6.9860 


-13.4700 


-0.0771 


-680.84003 


-681.55902 


sin 


0.20 


0.20 


-5.0490 


-14.1430 


-0.0574 


-687.51001 


-688.55798 


™Ge32 


2.50 


2.50 


-10.1990 


-6.7730 


-0.1831 


-610.54401 


-609.46198 


"Ge32 


1.00 


2.50 


-9.7040 


-7.4030 


-0.1692 


-617.96600 


-616.29401 


^^Geaa 


2.50 


2.50 


-9.4920 


-7.8830 


-0.1650 


-628.70801 


-628.55902 


^'Ge32 


1.00 


2.50 


-8.8660 


-8.6950 


-0.1935 


-635.50000 


-635.09198 


^*Ge32 


2.00 


2.00 


-8.6880 


-8.9680 


-0.1749 


-645.69202 


-644.26599 


'''Ge32 


1.00 


2.00 


-8.8110 


-9.3880 


0.1881 


-652.20398 


-650.90198 


^«Ge32 


1.50 


2.00 


-8.4500 


-9.9400 


0.1721 


-661.62598 


-660.47302 


^'Ge32 


1.00 


1.50 


-8.2590 


-10.3350 


0.1793 


-667.69800 


-667.63300 


^«Ge32 


0.75 


1.00 


-7.9020 


-10.8030 


0.1760 


-676.41998 


-675.31500 


™Ge32 


0.50 


0.25 


-7.6670 


-11.4350 


0.1660 


-682.12201 


-682.85199 


«°Ge32 


0.25 


0.25 


-6.7390 


-12.1790 


0.1530 


-690.15399 


-690.39801 


*iGe32 


0.25 


0.25 


-6.3290 


-12.7720 


0.1149 


-695.07599 


-696.47198 
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Table 7. Results of RMF calculations for various isotopes 
of Ni using the parameter set NL3. The values of the fixed 
proton and neutron pair-gaps (A„ and Ap) used in these 
calculations are identical to the values used in calculations 
with parameter set NL-SH. The remaining columns give 
the computed values of /i„, jip, and the differences between 
the calculated and experimental binding energies in table 
^ in MeV. Negative values mean that the calculated bind- 
ing energy is larger then the experimental binding energy. 



Isotope 



fj.p 



-BEexp. 



-B-Bnl-sh 
— BEexp. 



^«Ni28 


-11.044 


-5.692 


0.74099 


0.56998 


''Ni28 


-10.085 


-6.282 


0.60501 


0.45209 


«°Ni28 


-10.041 


-6.861 


0.48603 


-1.23132 


«^Ni28 


-9.527 


-7.556 


0.44602 


-1.46258 


«'Ni28 
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Fig. 2. The chemical potential of our RMF-calculations 
with the parameter set NL-SH ^inlRMF is compared with 
that of the BCK-model fJ.n\BCK and with the functional 
forms of the liquid drop model (LDM) given in Eq. (Q) and 
the finite range liquid drop model (FRLDM) ^IuIfrldm 
discussed in Sect. 4.2 for Mn isotopes 



Fig. 3. Experimental vs. calculated values of binding en- 
ergies using the parameter set NL-SH.. The calculations 
were carried out with the the pairing gap parameters listed 
in the Tables |, |, | and |. 

Fig. 4. Neutron chemical potentials of Ni isotopes calcu- 
lated within RMF theory based on the parameter set NL- 
SH and within the BCK-model are compared with the 
functional form of /i„ in Eq. (|l6| ) (solid line) and Eq. ^ 
(dashed line). 

Fig. 5. The chemical potentials ij-uIrmf (parameter set 
NL-SH) and ^uIbck for Ni. The values ^inlscK have been 
generated from the BCK-EOS using Sy = 30.34 MeV but 
with the same incompressibility K — 180 and ay — 16.0. 



Table 8. Comparison of the nucleonic chemical potentials 
for RMF and El Eid & Hillebrandt EOS in stellar collapse 
using data for a star of initial mass 15Mq. 



Model pio Ve Mean Nucleus ^„ (MeV) 



Fig. 6. Difference of proton and neutron chemical po- 
tentials fl\RMF (parameter set NL-SH) and filscK for Ni. 
The BCK- values have been generated from the BCK-EOS 
using Sy = 30.34 MeV, but with the same incompressibil- 
Mp (MeV)ity K = 180 and a„ = 16.0. 



RMF 0.37 
EH 0.37 



0.42 
0.42 



52 



-6.93 

-7.08 



-11.7 
-11.5 



Fig. 1. (a) Cumulative neutrino flux up to pio = 24.16 
g/cm^, computed with FuUer/BBAL EOS for M = 25 
Mq, D = 1 kpc and |A/gtP = 2.5 and later 0.1. (b) The 
spectrum in (a) folded with the detection cross-section for 
c.c. reaction Ve{d,pp)e~ in SNO. (c) The spectrum in (a) 
folded with the detection cross-section for v^ — e~ scat- 
tering in Super Kamioka. The projected lower limit E,^ is 
> 5 MeV for the S-K and SNO detectors. 



-if'ig. 7. The chemical potentials fJ.n\RMF and /inlscif, 
along with the functional forms of the liquid drop model 
(LD M) given in Eq. (^ and HuIfrldm discussed in 
Sect. 4.2 for Ni isotopes. This figure is especially signifi- 
cant because all the Ni isotopes are spherical nuclei. 



Fig. 8. The chemical potential of neutrons obtained with 
RMF theory based on the parameter set NL-SH is com- 
pared with the functional form of the Finite Range Droplet 
Model (FRDM) for Ni isotopes. 

Fig. 9. The chemical potentials HuIrmf and ^Jin\BCK, 
along with the functional forms of the liquid drop model 
(LD M) given in Eq. (|j) and iIuIfrldm discussed in Sect. 
\a for Fe isotopes. 



Fig. 10. The chemical potentials /JnlflA/ f an d ^inlscK, 

Eq. (|^) and iIuIfrldm discussed in Sect. 

topes. 



4.2 



for Ga iso- 
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Fig. 11. The nuclear binding energies using the spherical, 
incompressible LDM, the FRLDM (using volume, surface, 
coulomb and coulomb exchange terms only) and the ex- 
perimental values for Ni isotopes. 

Fig. 12. The nuclear binding energies using the spherical, 
incompressible LDM, the FRLDM (using volume, surface, 
coulomb and coulomb exchange terms only) and the ex- 
perimental values for Ge isotopes. 

Fig. 13. The difference of neutron and proton chemical 
potentials (/t — [in — /ip) for Ni isotopes, taking into ac- 
count the lattice correction, at a density of 10^^ g/cm'^. 
The figure also shows for comparison, the values of /i|flA/F 
without lattice correction, and the (lattice independent) 
LDM value for /< from Eq. (|). 
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